Observing Lense-Thirring Precession in Tidal Disruption Flares 
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When a star is tidally disrupted by a supermassive black hole (SMBH), the streams of liberated 
gas form an accretion disk after their return to pericenter. We demonstrate that Lense-Thirring 
precession in the spacetime around a rotating SMBH can produce significant time evolution of the 
disk angular momentum vector, due to both the periodic precession of the disk and the nonperiodic, 
differential precession of the bound debris streams. Jet precession and periodic modulation of disk 
luminosity are possible consequences. The persistence of the jetted X-ray emission in the Swift 
J164449. 3+573451 fiare suggests that the jet axis was aligned with the spin axis of the SMBH 
during this event. 

PACS numbers: 98.62.Js, 98.62.Mw, 98.62.Nx 



Introduction. The tidal disruption of a star by a su- 
permassive black hole offers a unique opportunity to 
probe the nuclei of otherwise quiescent galaxies. How- 
ever, the small number of candidate tidal disruption 
events (TDEs) makes it difficult to resolve theoretical un- 
certainties concerning their rates [iHil, super-Eddington 
accretion phase , and the period during which dissi- 



pation in shocks allows an accretion disk to form |10|, . 



An additional outstanding question about TDEs is 
whether or not they produce jets, as observed in many 
other accreting black hole systems. The past year has 
seen both the first theoretical models for TDE-associated 




jets [12, ll3| and the discovery by the Swift satellite of an 
intense, transient gamma- and X-ray flare from a galac- 
tic nucleus at z w 0.35 [IJL This flare has been ex- 
plained by multiple authors |15l4l7| as jet emission from 
a TDE aligned with our line of sight (although alternate 
hypotheses exist [S])- A second possible TDE-associated 
jet was also recently observed ^19.] . 

If such jet emission is common, then TDEs provide 
a unique probe of the physics of accretion and jet pro- 
duction in the vicinity of distant black holes' horizons. 
Specifically, it is unknown at present whether jets will 
align with the black hole spin vector, the disk angular 
momentum vector, or some other component of the mag- 
netic field geometry [23]. In most black hole accretion 
environments these directions are parallel, but the tran- 
sient disk of a TDE will generally have some tilt with 
respect to the SMBH equatorial plane. In this Letter 
we demonstrate that if jets from tilted TDE accretion 
disks align with the disk normal vector, they will gener- 
ally be expected to precess, often by observable amounts. 
Even absent the existence of a disk-aligned jet, or any jet 
at all, general relativistic (GR) effects will precess TDE 
disks with potentially observable consequences. 

Spin evolution of a tilted disk. Stars of mass A/* and 
radius i?* that pass within a radius 



FIG. 1. Geometry of the tidal disruption of a star by a spin- 
ning SMBH. Following disruption of the star near its pericen- 
ter passage, an accretion disk will form in the star's orbital 
plane. As the disk precesses, the angle /3 between the SMBH 
spin vector Jbh and the disk angular momentum vector Ldisk 
stays constant, but an associated jet may move relative to the 
observer's line of sight fobs- 



of a black hole of mass A/bh will be tidally disrupted, 
with half their mass immediately unbound from the black 
hole 21]. For black holes of mass Mbh > 10^ Af©, the 
tidal radius i?t is inside the Schwarzschild radius i?s and 
stars are swallowed whole rather than disrupted. The 
bound debris rapidly expands and cools so that its pres- 
sure is negligible and the approximation of geodesic mo- 
tion is accurate [TO]- The most tightly bound debris 
stream of a star disrupted at radius i?p returns in a time 



tfall « 50 M( 
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i?t =i?*(^^BH/A/*)l/3 
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where Afg = AfBH/lO^Afg, = R^/Rq and = Rp/Rs 
[1], although ifaii depends on the stellar density profile 
and can be evaluated more precisely by numerical simu- 
lations [ll]. After a small multiple of this time, stream- 
stream collisions circularize the returning gas and allow 
an accretion disk to form. In general, this transient ac- 
cretion disk will not lie in the black hole equatorial plane. 

An accretion disk inclined out of the equatorial plane 
of a spinning black hole by an angle (3 (assumed to equal 
the inclination of the stellar orbit before disruption, /3* 
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- see Fig. [T]) will be subject to Lense-Thirring torques 
with a strong radial dependence. For a thin disk |2S 



it is expected that the Bardeen-Petterson effect [2J, [2£ 
will induce a warp in the disk structure. However, for 
the thicker disks expected in many TDEs [1, [llj , simu- 
lations combining GR and magnetohydrodynamic effects 
(GRMHD) have shown that the disk precesses as a solid 
body rotator [26|, l27| . Such an accretion disk will pre- 
cess with a period Tp^ec = 27rsin/3 (J/r), where J is total 
angular momentum and t is integrated torque. A no- 
table feature of this formula is that Tpi-ec is independent 
of many disk model parameters, and depends only on the 
dimensionless radial surface density profile. 

The simulations mentioned above considered disks 
with a roughly constant surface density. For surface 
densities of the form S = , the precession 

timescale is 1261 



87rGMBH(l + 20 rV^-^r\'^+^{l - {n/r^f'^-^) 



c3(5-2C) 



a(l - (n/ro)i/2+<) 



(3) 

Here the disk inner and outer {Rq) edges have been 
normalized to units of Schwarzschild radii (n = Ri/Rs, 
To = Ro/Rs)- The variable a is the dimensionless black 
hole spin parameter, with values between and 1. 

Whether or not the disks associated with tidal disrup- 
tion flares approximately follow a surface density profile 
S = Ei(i?/fi!i)~^ is unclear. Ref. [8] presented a slim 
disk model for TDK accretion flows, with height H given 
by: 



H 
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Edd 
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where the function K is defined as 




Here / = 1 - {R-jRf/'^. M/MEdd is the ratio of the 
mass accretion rate to the Eddington rate assuming 10% 
accretion efficiency. 

However, this model is not suitable for use in calcu- 
lating Tpi-cc, as the zero-torque boundary condition used 
to calculate / leads to an unphysical singularity in S at 
i?i (S oc R^/'^K'^ / f). A different, numerical model was 
recently presented in Ref. (28.], in which axisymmetric 
disk equations were evolved with a time-dependent rate 
of mass input at the pericenter of disruption. This model 
led to a shallow decline of S with decreasing r after the 
arrival of the inner edge of the accretion flow at the inner- 
most stable circular orbit. Motivated by Refs. [s, 28|, we 



consider C, = —3/2,0, 1 in this paper. Tprec increases by 
a factor k, 7 when going from the C = 1 to the C, = —3/2 
model. 

The framework we followed is based on two under- 
lying assumptions: (i) & coherent accretion flow exists; 
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FIG. 2. Timescales for avoiding Bardeen-Peterson warping 
ithin (blue, top panel) and for establishing an accretion disk 
tcirc (green, bottom panel) as functions of the black hole mass 
Mbh- Dotted lines correspond to stars with a mass of 2Mq, 
solid lines to IMq and dashed lines to 0.5Mq (with a stellar 
mass-radius relationship adopted from Ref. [231 , p. 208). We 



take rio 



3 and Rp = 0.5-Rt, and conservatively plot tt 



for the outer edge of the disk, assuming Ro = 2Rp 



and (ii) the flow is not susceptible to Bardeen-Petterson 
warps {H / R > a, where a is the dimensionless disk vis- 
cosity parameter). Assumption (i) is only valid after a 
time tcirc ~ 'T-orbifaii, whcrc TT-orb is the number of or- 
bits required to circularize the most tightly bound debris 
streams fll|. A value of norb ~ 1-10 is often assumed 
in the TDE literature, but this quantity is poorly con- 
strained and could be higher for large a and /?*, where 
Lense-Thirring precession can delay the stream-stream 
collisions necessary for disk formation At later 

times, assumption (ii) will break down, as M/M-Edd de- 
clines and the disk becomes geometrically thinner. 

Adopting Eq. ^ for convenience, H/ R will fall below 
a after a time 
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/ Rs\ 
X-i R J 



where the function X ^ a and is X ^ a/{l — 8a^/3/). 
Also note that X^i — X/0.1. TDE disks will precess as 
solid body rotators during the time range tdrc < t < tthin 
as illustrated in Fig. [21 which shows that for Mbh 
IO^Mq (and any realistic Rp), solid body precession will 
occur for < 1 yr. 
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Angular momenta of returning debris streams. The 
evolution of the debris streams prior to their first return 
to pericenter has been studied in detail by Ref. [loj . 
The orbits of these streams, if non-equatorial, lack a con- 
stant orbital plane due to Lense-Thirring torques. The 
accretion disk is therefore fed by a supply of new gas 
with time-dependent angular momentum, which in turn 
evolves the direction of Ldisk- In. contrast to direct preces- 
sion of the accretion disk, we call this effect "differential 
stream precession," or DSP. Although we will compute 
numerical general relativistic solutions for the DSP, we 
can gain valuable intution from a simpler, lowest order 
estimate in the post-Newtonian limit. 

The angle by which the angular momentum vector of 
a debris stream will precess during an orbit of period T 
will be (/)orb(T) w Ar2sin(/3) = 27rsin(/3)(r/tLT), where 
AD, is the nodal precession and 



T 
2a 



c^Ajl - e^) 
GMbh 



3/2 



(7) 



is the Lense-Thirring precession period jSOf for a gas 
stream of semimajor axis A and eccentricity e. Defin- 
ing A0orb = (f>oihiT) — (/)orb(oo) as a mcasurc of the DSP, 

A(/.orb = 47rasin(/3)(2rp)3/2((i + e)-^/^ _ 2-^/2), (8) 

which Taylor expands in the late-time, Rp/A <§; 1 limit 

to A0orb « 1.7 sin(;3)arp ^/V,M6-i(tAfaii)-'/^. 

Although Eq. ([8]) is not exact, it provides a valuable 
insight: the DSP is largest for low-mass, rapidly spinning 
SMBHs that disrupt stars with deeply plunging, inclined 
initial orbits. At early times the disk viscous timescale 
ivisc ^ ifaii [1] so Eq. ([U approximates the angular evo- 
lution of I/disk- We do not expect A0orb > 1° after the 
establishment of a steady accretion flow {t > Stfaii) for 
any TDEs with solar- type stars and Mq > 2, although 
these constraints relax for stars with r, > 1. 

To obtain an exact solution for the time evolution of 
angular momentum in the returning debris streams, a 
GR calculation is needed. We numerically integrate the 
Kerr geodesic equations following the formalism of Ref. 
3l|. We assume a flat distribution of debris mass with 
specific Newtonian energy a spread in that energy of 
3GMbhR*/RI Hi, and obtain constants of integration 
for each debris stream by transforming the initial condi- 
tions {E, Rp, /3} to {Eqr, iz, Q} {Egr, iz, Q are specific 
energy, z-component angular momentum, and Carter's 
constant for Kerr metric test particles) . Good agreement 
with Eq. ([U is shown in Fig. |3l 

Observational implications. We have shown that the 
Lense-Thirring effect will cause the direction of a TDE 
disk's angular momentum vector to evolve in time. Di- 
rect precession of the accretion disk is the dominant ef- 
fect, but in some cases DSP can cause a significant non- 
periodic evolution in Ldisk- The precession of the ac- 
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FIG. 3. The angular shift A(/>orb- The thick curves illustrate 
the disruption of a solar- type star with Mbh = 10^ M©, a = 
0.8, and rp = 13; the thin curves are the same but with rp = 3. 
The blue dotted lines are Eq. ((8]), while the green solid lines 
are numerical geodesic solutions. The curves do not extend 
prior to i = tfaii, and are normalized by sin/3. 



cretion disk will modulate the observed disk luminos- 
ity at least by a factor of cos('!/'), and lead to periodic 
pulsations of the associated transient quasar. This peri- 
odic modulation could in principle be extracted from the 
Fourier decomposition of a TDE lightcurve, but perhaps 
a more promising avenue for detection lies in the fraction 
of events for which the disks will precess into an edge-on 
phase. This could reduce the observed disk flux by ~ 2 
orders of magnitude while simultaneously reddening the 
peak emission frequency fill . Even in the absence of jet 
emission, observations of a "blinking" TDE flare could 
provide strong evidence of precession and allow both a 
and the disruption parameters to be constrained. 

The most exciting possible consequence, however, is 
precession of jets associated with TDE disks. If we as- 
sume that relativistic jets in tilted accretion systems align 
with Ldisk, narrow jets will precess out of the observer's 
line of sight in a small fraction of Tprcc- Continuous ob- 
servation of a jet for a relatively short period of time, iobs, 
allows very strong constraints to be placed on combina- 
tions of a and disruption parameters such as rp and /3* 
(assuming still that (3 — /3^,). Alternatively, repeated ob- 
servation of TDE-associated jets could serve as evidence 
that jets align with Jbh or an aspect of the magnetic 
field geometry, provided that sufficient non-precession is 
observed. We note that the DSP, though generally sub- 
dominant, can in some cases cause very rapid precession 
(up to ~ 0.1°/min) at the onset of the flare (Fig. [3]). 
If an associated jet is aligned with Ldisk, this will lead 
to a brief, nonrepeating transient which could fake an 
unusually long gamma ray burst provided 0jct ^ 1°. 

To provide a concrete example of the above consider- 
ations, we consider the tidal disruption candidate Swift 
J164449.3-H573451, for which Ref. [H inferred the fol- 
lowing relevant disruption parameters: Mbh ^ 10"" — 
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FIG. 4. Regions of a-/? parameter space that can be excluded 
by continuous observations of a TDE jet with the inferred 
parameters in Ref. [ISj] and = 0. The sohd curves show 
contours of constant fobs = Tprec x 2(Sjct/10~^'^)/(27rsin/3): 
the maximum number of days it would take for a jet initially 
in the observers' line of sight to precess off-axis, with the jet 
opening angle normalized to 10~^'^. We take Ro = 2Rp and 
Ri = 3Rs- Regions of parameter space to the right of the 
thick red contours can be excluded for the Swift TDE jet, 
which exhibited bright X-ray emission for over two weeks. 
The 14 day contours for ( — —3/2 and C = 1 are shown with 
black dotted and dashed lines, respectively. The effect of the 
DSP is negligible for these parameters, and neglected here. 



to the observed rate of jets over the period of the Svifift 
mission, and the Eddington hmit of the SMBH) . 

Figure 4 shovifs the resulting constraints on the joint a- 
(3 parameter space of this TDE if wc take Mg = 1 . Since 
the bright X-ray emission from Swift J164449. 3+573451 
persisted for over two weeks, at least one of the follow- 
ing statements must be true: (i) the value of a is ex- 
tremely low, < 10-2 (10-1 if C = -3/2); (ii) the initial 
orbit of the disrupted star was tightly aligned to within 
~ 0jet of the black hole equatorial plane; or (in) the jet 
emission was not aligned with the disk spin axis. The 
first possibility would represent an unusually low value of 
black hole spin and could be excluded if the Blandford- 
Znajek mechanism was responsible for jet launching (32i] . 
while (ii) requires that there will be a larger abundance 
of somewhat shorter events. Since such flares are not fre- 
quently observed, the persistent X-ray emission in Swift 
J164449. 3-1-573451 suggests that its jet was aligned with 
the steady spin axis of the black hole rather than with 
its precessing disk. Future GRMHD simulations can test 
this inference from first principles. The detection of ad- 



ditional TDE-associated jets in future surveys would test 
the statistical robustness of this conclusion. 
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